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FACTORIZATION OF BANDED PERMUTATIONS 



GRETA PANOVA 



CN Abstract. We consider the factorization of permutations into bandwidth 1 

T" H permutations, which are products of mutually nonadjacent simple transposi- 

Q_^ tions. We exhibit an upper bound on the minimal number of such factors and 

CN thus prove a conjecture of Gilbert Strang: a banded permutation of band- 

width w can be represented as the product of at most 2ui — 1 permutations of 
bandwidth 1. An analogous result holds also for infinite and cyclically banded 
permutations. 
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1. Introduction 



Computational efficiency very often requires us to represent matrices as products 

of certain special, easily computable, matrices using as few factors as possible. 

Matrices of bounded bandwidth are often seen in practical applications. In p\ and 

[6] Gilbert Strang shows that when a matrix and its inverse are of bandwidth w, 

it can always be represented as a product of 0(w 2 ) such matrices of bandwidth 

£\| w = 1. In particular this bound is independent of the size of the given matrix. He 

J> also conjectures that for permutation matrices this bound is actually 2w — 1. In this 

paper we will prove this conjecture. This conjecture has been proven independently 

^^ later also by Albert, Li, Strang and Yu in [lj and Ezerman and Samson in 13]. 

A matrix of bandwidth w is a matrix A, whose nonzero entries lie within 

distance w from the main diagonal: A\j = whenever | i — j | > w. In particular, 

a banded permutation matrix P is a — 1 matrix with exactly one 1 in each row 

and column and such that Pjj = if | % — j \ > w. The matrix P corresponds to 

,-H the permutation tt defined as 7r^ — j for Pij — 1 and vice versa. So tt is of width w 

lj if 1 7Tj — i | < w for every i. 

• ^h Our main result is the following. 

>< 
$_j Theorem 1. For any permutation n, let M = {ttj — i\i < j,ni > ttj} and let 

m = #M be the number of different elements in the set M. Then there exist at 

most m bandwidth 1 permutations pP , i.e. | p? — i | < 1, such that it — p 1 p 2 

If 7r is a permutation of bandwidth w then for the elements of M we have that 
7Tj — i < iTi — i < w and ttj — i > ttj — j > —w, so M C {— w + 1, . . . ,w — 1}. Thus 
M has no more than 2w — 1 elements and the conjecture follows immediately. 

Corollary 2 (Strang's conjecture). If it is a permutation of bandwidth w then 
there exist at most 2w — 1 bandwidth 1 permutations whose product is it. More- 
over, the bound 2w — 1 is exact. 
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A possible extension of banded matrices, also considered by Strang in m\, are 
infinite and cyclically banded matrices. Cyclically banded matrices are n x n ma- 



trices A, such that Ai 



if w < 



< n 



Here any matrix would 



have width n/2 and so we will require that w < n/2. In Section^ we consider the 
analogous question referring to their factorization into cyclically banded matrices 
of bandwidth 1. 

We will use the notion of reduced decomposition of a permutation and its vi- 
sualization called a wiring diagram. We would like to thank Alex Postnikov for 
suggesting their use. The proofs will rely on the construction of special wiring 
diagram we call a hook wiring diagram. 

2. Hook wiring diagrams 

We will consider the simple generators of the symmetric group S n as a Cox- 
eter group. We will represent a permutation as a certain product of such simple 
transpositions which will be grouped into the desired bandwidth 1 factors. 

A simple transposition Si — (i, i + 1) exchanges the ith and (i + l)st element. As 
an element of the symmetric group S n , Si is equal to the permutation 1, 2, . . . , % — 
1, i + 1, i, i + 2, . . . , n. Note that i determines the transposition, we will say that i is 
the index of the transposition Si. A reduced decomposition of ir is a product 
Si±Si 2 ■ ■ ■ Si t — it of such transpositions of minimal possible length 1; see p] for the 
general facts. It follows by inspection of the possible cases that every bandwidth 1 
permutation is a product of mutually nonadjacent simple transpositions s i; where 
two transpositions are adjacent if their indices are consecutive numbers. 

A wiring diagram (originally appearing in El) of a reduced decomposition 
s% 1 8i 2 ■ • • s^ — ir is a planar configuration of n (pseudo-)lines L\, . . . , L n between 
two columns of the numbers 1, 2, . . . , n with the following properties: 

• Line Li starts at i and ends at m. 

• No two lines intersect more than once and no three lines intersect at a point. 
Each wiring diagram depicts a reduced decomposition Si ± Si 2 ■ ■ ■ Si t — ir via the 
following correspondence. Through every intersection between the lines draw a 
perpendicular "dashed" line (as in Figure [I]). Counting from left to right assign 
to the fcth dashed line the simple transposition Si k = (ik,ik+i), whose index ik is 
equal to 1 plus the total number of (pseudo-)lines L r that cross that dashed line 
above the intersection point. 




Figure 1. Wiring diagram of 7r = S2S4S3S1S4 = 25143. 

Figure [l] shows the wiring diagram for ir = (2,3)(4,5)(3,4)(1,2)(4,5). Notice 
that line Li "carries" the number i. A thin vertical slice of a wiring diagram 
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represents an intermediate permutation with the position of i being the relative 
vertical position of Li with respect to the other lines at this slice. Two lines 
crossing simply means that we exchange two adjacent entries. The number of lines 
vertically above that crossing plus 1 is exactly the index i of the corresponding 
simple transposition s, = (i,i + 1). 
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Figure 2. Hook diagram for the permutation 7r = 5,4,7,1,9,2,3, 
10,8,6 of bandwidth w = 4. The dashed lines depict the diagonals 
c — r = —3,-2,-1,0,1,2,3. The numbers at the intersections 
are the indices of the corresponding simple transposition, e.g. 2 
corresponds to S2 — (2,3). 



For any permutation n we can also draw (see Figure [2]) what we'll call a hook 
diagram. Consider a square grid bounded by (0, 0) in the top left corner and ver- 
tical and horizontal rays marked with 1,2,... going down and to the right following 
the indexing convention for matrices; so that a point of coordinates (r, c) is at the 
r-th row (counting from the top) and c-th column (counting left to right). Place a 
dot at the points (i, 7r») on the grid and connect (i, 7Tj) with (i, 0) and (0, m) by two 
segments. This way the dots would be at the places of the ones in the permutation 
matrix of n and each i will be connected to the corresponding 7r^ by a hook with 
corner at the dot (i,TTi). 

Notice that a hook diagram turns readily into a wiring diagram by extending 
the horizontal segments through (0, i) and the vertical segments through (J, 0), then 
rotating by —45° as shown in Figure [3j The line Li would be the rotated extended 
hook through the points (0, i), (iri, i), (wi, 0). To determine the index of the simple 
transposition corresponding to a crossing of Li and Lj we need to count the number 
of lines in a thin strip vertically above that crossing in the rotated extended diagram. 
Assume i < j, so since Li and Lj cross we must have 7r^ > ttj. The number of lines 
above the crossing and hence the index will be r — i + ftj ~#{p : P < i, ftp < ttj} + 1, 
so the transposition is s r = (r, r + 1). 

Proof of Theorem U\ Draw the hook diagram of 7r and interpret it as a wiring di- 
agram as described above. We can now read off a reduced decomposition from 
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Figure 3. The wiring diagram obtained from the hook diagram 
for the permutation 7r = 5, 4, 7, 1, 9, 2, 3, 10, 8, 6. The numbers at 
the intersections indicate the indices of the corresponding simple 
transposition. 

the hook-wiring diagram as follows. To every intersection of two hooks assign the 
number i of the corresponding simple transposition as explained above. 

Let a k be the product of the transpositions on the fc-th diagonal, i.e. c — r — k 
(see Figure p2j), where k = — n+1, . . . ,n— 1. We have a k = Si t ■ ■ ■ s,, where i\,...,ii 
are the numbers(indices) at the crossings on the fcth diagonal. These numbers are 
at least 2 apart, so we have a k — i if i £• {ii, i\ + 1, . . . , ii, ii + 1}, of. = ij + 1 
and a\ , 1 = ij for j = 1, . . . ,1. Then er fc is of bandwidth 1. The diagram gives a 
reduced decomposition of 7r, where the transpositions corresponding to intersections 
on the same diagonal c — r — k appear before the ones on the next diagonal, while 
the relative order of these transpositions on the same diagonal does not matter 
since they commute. Since their product is a k we have 7r = <j~ n+l a~ n+2 ■ ■ ■ a n ~ 1 . 
Notice that a k is not trivial if and only if there is an intersection on this diagonal. 
There is an intersection between Li and Lj, (i < j), if and only if 7r^ > ttj, and the 
intersection point is (i, ttj). Thus diagonal k has an intersection if and only if k G M 
and the number of nontrivial bandwidth 1 factors is equal to the cardinality of M, 
m = #M. If M = {ki, . . . , k m }, then the nontrivial factors are a ki for i = 1, . . . , m 
and we set p l = <j ki to obtain 7r = p x p 2 ... as desired. □ 

Proof of Corollary^ As remarked earlier, if it is banded then M C {— w+1, . . . , w— 
1}, so #M < 2w — 1 and Strang's conjecture follows from Theoremfll 

To show that 2w — 1 is the exact bound, consider the permutation a — (w + 
l)(u) + 2) . . . (2w)123 ... w ... of width w, where the last . . . mean the identity o~i = i 
for i > 2w. Before we show that this particular a cannot be factored into less than 
2u> — 1 permutations of bandwidth 1, we need to make a few general observations. 

For any permutation 7r, let k be the minimal number for which 7r = 7P 1 ) • • • ir( k \ 
where 7rW are permutations of bandwidth 1. Then there exists a reduced decom- 
position of 7r = sljslj ■ ■ ■ s\ , such that n^ is the product of the ith block of 
transpositions. We will show that there is always such reduced decomposition by 
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8 8 

9 9 

Figure 4. Any wiring diagram of a = (w + 1) . . . (2w)12 .. .w ... 
is homotopy equivalent to this one. Here w — 3 and a = 
456123789. 

decreasing the number of simple transpositions in it. Writing tt^' — Si ± ■ ■ ■ Si m as 
a product of transpositions we still have a decomposition of it into simple trans- 
positions. We can depict this decomposition graphically like a wiring diagram, 
without requiring that two lines intersect at most once. The assertion that the 
decomposition of 7i"W = s^ ■ ■ ■ s im is not reduced is equivalent to two lines L' and 
L" intersecting at least twice at places r and p corresponding to Si r and si . Let 
L' = A'B'C" and L" = A"B"C" where A, B, C are the portions of the lines ob- 
tained after cutting at the two intersections. Substituting L' and L" with A'B"C 
and A" B'C" respectively gives us another wiring diagram of 7r for the decomposi- 
tion 7r = Si ± ■ ■ ■ Si r ■ ■ ■ si ■ ■ ■ Si m . Removing Sj r and Si from the 7p''s to which they 
belonged gives another factorization of it into at most k permutations of width 1 
with a smaller number of simple transpositions. Continuing this way we will reach 
the length / of n forcing the underlying decomposition into simple transpositions 
to be reduced. 

We can thus assume that our particular a = a^ 1 ' ■ ■ ■ cr* > gives a reduced decom- 
position. Consider its wiring diagram as depicted in Figure [4] since the transposi- 
tions in each a^' are nonadjacent we can draw the corresponding intersections on 
the same vertical line. Thus every path from some i to some Oj will pass through 
at most k intersections. 

Notice that any wiring diagram of a can be deformed (is ambicntly isotopic) to 
a w x w grid rotated 45° like the diagram on Figure [4] Then every path joining 
w + 1 with a\ = w + 1 has 2w — 1 intersection points and so k > 2w — 1. □ 

Since our proof is constructive, it leads to an algorithm for the decomposition: 
find the intersection points in the hook diagram and group them according to the 
diagonal to which they belong. 

Let Ik be the set of intersection points on the fcth diagonal. Assume the inverse 
permutation 7r _1 is known. Then the procedure is as follows: 

For i from 1 to n: 

P ■= TTj 

For j from 1 to p — 1 : 

If ttJ 1 > i, then Ij_j <— (i,j) 
In order to determine which transposition these intersections correspond to, no- 
tice that the number of lines L r intersecting the segment between (i,j) and the 
origin, and thus the index of the transposition, is i—l+j — 1 — #{i | t < i , ir t < j}, 
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which we can count within this algorithm also. Let s[i,j] = #{£ \ t < i , ir t < j}, 
then 

s[i,j + l] = s[i,j] + ((irr 1 <i)), 

where ((statement)) denotes the logical value 0/1 of the statement. 

3. Infinite and cyclically banded permutations. 

We now consider an extension of banded matrices and their corresponding per- 
mutations n = 7i"i . . . ir n , as defined in Section [T] In this case bandwidth 1 encom- 
passes more permutations and thus allowed factors. We have a simple transposition 
so = s n exchanging -k\ and 7r„ , and corresponding to the bandwidth 1 cyclic matrix 
A with A\ n = A n i = 1. The other additional factor is the shift S, acting by cyclic 
shift on 7r as Sir = 7r 2 7T3 . . . 7r„7ri, with corresponding matrix S given by Sa+i = 1 
and otherwise. A bandwidth 1 permutation is either a shift S or S" 1 or a product 
of pairwise nonadjacent modulo n simple transpositions sq, • • ■ , s n _i. 

...0 100000 0... 
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01 000000 00 

10 000000 00 

00 010000 00 

00 000100 00 

00 001000 00 

00 000000 10 



00 000001 00 
...0 000010 0. .. 

Figure 5 . The matrix A of the cyclic banded permutation n = 
652431 of bandwidth w = 3 and the corresponding infinite B = 
4(A). 



A cyclic banded matrix can also be interpreted as a doubly infinite periodic 
matrix of period (n, n) in the following way. For any nx n matrix A define 4>(A) = 
B, where B is a doubly infinite matrix given by B^j = Aj ( mo( j n \j ( mo( j n \ for 
I i — j | < fi/2 and otherwise, see Figure [5] for an example. The map <f> from cyclic 
banded matrices of bandwidth w to banded doubly infinite matrices of period (n, n) 
is an isomorphism. We will thus consider the problem of factorization of banded 
doubly infinite periodic matrices and their corresponding infinite permutations. 

We can now form the hook wiring diagram of the infinite periodic banded per- 
mutation matrix as in the previous Section and as it is shown on Figure [6J Again 
the intersections on each diagonal correspond to a bandwidth 1 factor and we can 
proceed to "untangle" them, i.e. multiply by simple transpositions. Let a 1 be 
the product of the simple transpositions corresponding to intersections on the «th 
diagonal, where each transposition is represented by its infinite periodic matrix 
and thus the number of these transpositions will be finite. After multiplication by 
a 1 ■ ■ ■ a 2w ~ 1 we will have a trivial wiring diagram where no lines intersect. Unlike in 
the finite case where a monotone bijective map [1, . . . , n] — > [1, . . . , n] must be the 
identity, in the infinite case of Z these could be shifts, so the permutation/matrix 
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Figure 6. The hook diagram for 7r = 652431, produced from its 
infinite matrix B, the corresponding wiring diagram which gives 
factors a 1 = si,er 2 = sos 4 ,cr 3 = Sj and the resulting trivial wiring 
diagram after "untangling" giving a shift S 1 . 



corresponding to a trivial wiring diagram will be S k for some k. We will now show 
that k < w. 

Define the shifting index of an infinite periodic permutation matrix B as fol- 
lows. Let p := #{Bij = 1 1 i < n and n + 1 < j} and q :— #{Bij = 1 1 n + 1 < 
i and j < n}. Alternatively these are the the number of ones in the upper right 
wxw triangle of the original A and the number of ones in the lower left such triangle. 
Let the shifting index of B be si(_B) := p — q. We have that si(id) = 0, si(S) = 1, 
si(5" 1 ) = -1. We also have that si(SB) = 1 + si(-B) and si(S , - 1 B) = -1 + si(B), 
since S acts on B by shifting one row upwards and the entry 1 on row n + 1 mov- 
ing to row n either decreases q by 1 or increases p by 1. Thus every permutation 
matrix factorizes uniquely as B = S S1 ^B, where si(B) = 0. Moreover, for any 
simple transposition Si and its corresponding infinite matrix Ei, we have si(.E;) = 
and si(EiB) = si(_B). This can be checked by inspection, i = (mod n) is the only 
nontrivial case, where it is still obvious that p — q is preserved after switching row 
n and n-\- 1 in B. Also SEiS^ 1 = Ei-\, where the indexing is modulo n again. 

Using the usual inversion index we can show by induction on it that every such 
B is the product of shifts and simple transpositions (note that any permutation 
matrix without inversions is a diagonal of Is, i.e. S k for some k). Thus we can 
always write B = S si ( B) E n ■■■E il . This also shows that si(AB) = si(A) + si(B), 
under repeated application of EiS = SEi + \. 

For a matrix B of bandwidth w we have that p < w and q < w, thus — w < 
si(B) < w. Moreover, B = a 1 ■ ■ ■ a 2w ^S k , so si(B) = si(cr 1 ) + - • ■+si(a 2w - 1 ) + k = 
k, and thus \k\ < w. We thus have the following analogue of Strang's conjecture, 
Corollary [2j 

Theorem 3. Let n 6 S n be a cyclic banded permutation of bandwidth w, i.e. 
| tt^ — i | <!oor I ni — i I > ii — w, or alternatively an infinite periodic permutation 
whose matrix is doubly infinite periodic with only nonzero entries in | i — j \ < w. 
Then n — a 1 . . . <j 2w ~ l S k : , where a % is a cyclic permutation of bandwidth 1, product 
of nonadjacent simple transpositions, and S is the cyclic shift by 1 with \k\ < w. 
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